Trapped ion chain as a neural network 
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We demonstrate the possibility of realizing a neural network in a chain of trapped ions with 
induced long range interactions. Such models permit to store information distributed over the 
whole system. The storage capacity of such network, which depends on the phonon spectrum 
of the system, can be controlled by changing the external trapping potential and/or by applying 
longitudinal local magnetic fields. The system properties suggest the possibility of implementing 
robust distributed realizations of quantum logic. 



The Cirac-Zoller trapped ion computer has become 
one of the paradigmatic models to implement a quan- 
tum computer 0. Recently, spectacular experimental 
progress in both realization of simple algorithms and im- 
plementation of quantum logic has been achieved using 
a few ions qubits (c.f. @). 

Still, serious experimental challenges remain on the 
way to a large-scale ion trap quantum computer. One 
should stress that while the achievement of an all-purpose 
quantum computer in the next years seems to be diffi- 
cult, one can be quite optimistic about the applications of 
chains of trapped ions as quantum simulators. Recently, 
it has been shown that long range pairwise interaction 
between individual spins is induced in an ion trap when 
introducingan additional state-dependent force acting on 
the ions [j, |5| . The properties of such ion spin molecules 
can be designed by variation of the parameters charac- 
terizing the ion trap. Later, a state-dependent optical 
force has also been shown to evoke long range couplings 
and has been proposed to simulate general Heisenberg in- 
teractions and to investigate quantum phase transitions 
in ion traps . Here we show that ion spin systems can 
serve to realize a neural network (NN) model. Neural net- 
works are a prototype model of parallel distributed mem- 
ory 0,13, an d have been intensively studied by physicists 
since the famous paper by Hopfield . These disordered 
systems with long range interactions typically present a 
large number of metastable (free) energy minima, like in 
spin glasses 8]. These states can be used to store infor- 
mation distributed over the whole system. The patterns 
stored by the network have large basins of attraction in 
the thermodynamical sense which means that even fuzzy 
patterns are recognized as perfect ones. For this reason, 
attractive NN's can be used as associative memory. At 
the same time, NN are robust, so that destroying even a 
large part of the network does not necessarily diminish 
the network performance. The above listed properties 
make NN's interesting for distributed quantum informa- 
tion storing and processing where a quantum bit would 
correspond to two collective modes of the ion's chain in- 
stead of two internal states of each individual ion. Differ- 



ent approaches to exploit the potential use of NN models 
for quantum information processing have been already 
discussed in the literature [1(1 ll l| . In we have shown 
that it is possible to create next neighbour entanglement 
in an ideal NN model. Here we propose, for the first time, 
a feasible implementation of NN using a chain of trapped 
ions. 

Before explaining the details of NN implementation we 
first remind the readers the main features of the Hopfield 
model ISJ, and discuss the effective Hamiltonian derived 
m The similarities between both models sug- 

gest the possibility of using a chain of trapped ions as 
a NN. To demonstrate this assumption we find the stor- 
age capacity of the system and its robustness which are 
the most appealing features of NN for distributed quan- 
tum information. Thus, the question of ergodicity and, 
therefore, the ability of the system to perform an asso- 
ciative memory is not considered here. We show that 
the storage capacity, which is determined by the phonon 
spectrum of the system, can be controlled by modifying 
the shape of the external axial trapping potential and/or 
by applying longitudinal magnetic fields. Although the 
storage capacity is a classical property of the network, 
spin ion systems permit also to study quantum neural 
networks by using a transverse magnetic field or an op- 
tical field that effectively simulates an additional trans- 
verse magnetic field. Transverse magnetic fields should 
permit tunneling processes between stored patterns and 
to realize for example of quantum stimulated annealing 
0, . This study is beyond the scope of this letter and 
will be treated elsewhere. 

Following the models of Hopfield and Little 01 > a 
neuron can be viewed as an Ising spin with two possible 
states: an "up" position (S = +1) and a "down" position 
(S = —1) depending on whether the neuron has or has 
not fired an electromagnetic signal, in a given interval of 
time |3- The state of the network of N neurons at a 
certain time is defined by the instantaneous configuration 
of all the spin variables {Si} at this time. The dynamic 
evolution of these states is determined by the interactions 
among neurons, Jh, which are symmetric, so, for any pair 
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of neurons, Jy = Jji. Also full connectivity is assumed, 
that is, every neuron can receive an input from any other 
one, and can send an output to it. The Hamiltonian of 
the system can be written as: 
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where h corresponds to an external magnetic field. The 
interactions are determined by the patterns or configura- 
tions of spins to be stored in the network. These patterns 
will be learned if the system is able to accommodate them 
as attractors, implying that a large set of initial config- 
urations of the network will be driven dynamically (af- 
ter sufficient enough time) to those patterns. A possible 
choice of the interactions is: 
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with i ^ j. The p sets of } = ±1 are the patterns that 
we wish to store. The network will have the capacity of 
storage and retrieval of information if the dynamical sta- 
ble configurations (local minima) reached by the system 
{Si} are correlated with the learned ones Despite 
the fact that the interactions have been constructed to 
guarantee that certain specified patterns be fixed points 
of the dynamics, the non-linearity of the dynamical pro- 
cess may induce additional attractors, the so called spu- 
rious states. 

Recently it has been shown that the Hamiltonian de- 
scribing a linear chain of harmonically trapped ions ex- 
posed to a magnetic field gradient Q or interacting with 
convenient laser fields [6( can be transformed into an ef- 
fective spin-spin Hamiltonian with long range interac- 
tions mediated by the collective motion of the ions: 
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where: 
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with a = x,y and z being the spatial directions, i,j la- 
beling the ions, a being the Pauli matrices, FX the forces 
in the a direction, m the mass of the ion, ui a ,n the angu- 
lar frequency of vibrational mode n. Mf n are the unitary 
matrices that diagonalize the vibrational Hamiltonian: 
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where nfj are the elastic constants of the chain ^|] . The 
coefficient M" n gives the scaled amplitude of the local 
oscillations of ion i around its equilibrium position when 
the collective vibrational mode n is excited. Thus, the 



eigenvectors of M describe each ion's contribution to a 
given vibrational mode while the eigenvalues provide the 
oscillation frequencies, u> a .n of the collective modes. 

The external trapping frequencies are chosen such that 
the laser cooled ions crystallize in a linear chain (i.e., 
Vl x = fly ^S> VI g) and the external forces are such that 
on the z— axis F£ — Ff 4 — holds, and the index a 
is dropped. From now on we will consider the case of 
a chain of ions interacting with external laser fields 
which allow us to consider site independent forces and 
magnetic field, F %j = F, B[ = B' 0. Thus, Eqs.lQJ 
and @ have the same form and the possibility to im- 
plement a classical neural network using a linear chain of 
ions arises. We substitute the Pauli matrices in Eq.@ by 
Ising spins S = ±1, where the effective spin corresponds 
to the internal state of the ion. Nevertheless, there are 
some important differences between both models. First, 
in the Hopfield model, the interactions given by Eq. (J2J) 
are determined by the patterns to be stored {£f } = ±1 
while in the trapped ion chain, the interactions are fixed 
by the coefficients M^ n that do not necessarily equal ±1. 
Second, in Eq.(J2J), p corresponds to the number of pat- 
terns to be stored, which in the limit of large number of 
spins N is bounded from above by p — 0.14/V 0- In 
our model, the sum (Eq.(0J) extends to the total num- 
ber of vibrational modes which is larger than the total 
number of stored patterns. And finally, in the Hopfield 
model all the patterns have the same weight while in the 
ion chain each vibrational mode (Eq.ipjl) is weighted by 
1 Due to this fact, to reproduce as close as possible 
a NN model, the vibrational modes should be ideally al- 
most degenerate in energy. Moreover, the corresponding 
patterns must have large basins of attraction, i.e. they 
should correspond to sufficiently different configurations 
of spins so that each configuration is dynamically recov- 
ered even if several spins are randomly flipped. These 
requirements are not achievable with an harmonic trap- 
ping potential due to the fixed frequency ratio between 
the two lowest vibrational modes • The desired degen- 
eracy can be achieved either by engineering an adequate 
external magnetic field or by using a trapping potential 
of the form of | x | 7 with 7 < 0.8 (see Fig^l. 

To check the feasibility of implementing a NN model 
in these ion spin systems, we calculate first the phonon 
spectrum, impose the learning rule (Eq.Q), and probe 
the correlation between the dynamically stable configu- 
rations and the stored ones. Therefore, the first step con- 
sists in finding the equilibrium positions of the ions and 
calculate the vibrational modes and frequencies by using 
a standard diagonalization procedure. Once the com- 
plete vibrational spectrum of the system is known, we 
calculate the interactions among spins J^, and evaluate 
the energy of the system. The dynamical thermalization 
of the system is simulated using standard Monte Carlo 
techniques. 

We fix the initial spin configuration by mimicking the 
signs of a given vibrational mode. Starting with this 
spin configuration, we flip randomly r spins and let the 
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FIG. 1: Ratio between the frequencies of the second and the 
first vibrational modes as a function of the exponent of the 
trapping potential, 7 



system evolve towards the equilibrium situation assum- 
ing a noiseless scenario. If the system recovers the ini- 
tial configuration, the configuration is stable under the 
flip of r spins. The number of initial spin flips deter- 
mines the initial overlap, defined as rrij = (N — r)/N. 
After dynamical evolution, the final overlap is given by 
nij = (N — s) /N where s is the number of spins that dif- 
fer from the initial configuration. We repeat this process 
over M initial different random r spin flips to evaluate 
statistically the final overlap with the initial configura- 
tion: nif = ('}2^ = i'mjnj)/M, rij being the number of 
times that the system reaches the rrij configuration. The 
value of the initial overlap for which significant decrease 
of the final overlap occurs, is a good measure of the size 
of the basin of attraction of the corresponding pattern. 
For the harmonic trapping potential, the two lowest vi- 
brational modes are the center of mass (all spins parallel, 
with frequency u>i = fi z ) and the breathing mode (half 
up, half down, with L02 = V3fl z ) E3- We find that the 
center of mass mode is stable up to the flip of half of the 
spins, while the breathing mode is already unstable under 
the flip of a single spin. Therefore, the only configuration 
that can be recovered is the center of mass mode, corre- 
sponding to all spins aligned up or down. To increase the 
storage capacity of the network one can engineer the ap- 
propriate spatial dependence of a longitudinal magnetic 
field to stabilize one additional given pattern. We apply 
this technique to make the second pattern stable. The 
results are displayed in Fig [3 showing the basin of attrac- 
tion of the two first modes for a chain of (a) 20 Ca + ions 
with B = 7.7 x 10~ 31 Tm 2 C/s parallel to each spin of the 
second mode and (b) 40 ions for B = 6.5 x 10~ 31 Tm 2 C/s 
parallel to each spin of the second mode. The results are 
averaged over 500 configurations and with random spin 
flip and random update. By increasing (diminishing) the 
external magnetic field the stability of the first pattern 
decreases (increases), while the second pattern exhibits 
the complementary behavior. In Fig^a) we show that 
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FIG. 2: Basin of attraction of the two lowest energy modes 
of an harmonically trapped (Q z = 1MHz) chain of ions in 
the presence of an external longitudinal magnetic field, (a) 
20 Ca+ ions with B = 7.7 x 10" 31 Tm 2 C/s parallel to each 
spin of the second mode and (b) 40 Ca + ions with B = 6.5 x 
W~ sl Tm 2 C/s parallel to each spin of the second mode. Black 
squares (triangles) correspond to the lowest (second) mode. 



with 20 ions, it is possible to stabilize completely the first 
and second pattern simultaneously up to the spin flip of 
3 ions while the probability of recovering the initial pat- 
tern for four initial random spin flips is 99.2% for both 
patterns. Increasing the number of ions to 40, FigEfb), 
the basin of attraction increases up to the random flip 
of 9 ions with 100% of probability of recovering of the 
initial patterns, 99.7% (99.5%) of probability of recover- 
ing the first (second) pattern for the initial random flip 
of 10 spins and 98.5% (99%) of probability of recovering 
the first (second) pattern for the initial random flip of 11 
spins. 

The second approach used to increase the storage ca- 
pacity of the network consists in modifying the spatial 
dependence of the trapping potential. We consider an 
external potential of the form V(x) — A \ x | 7 with an 
even number of ions in order to avoid any possible sin- 
gularity of the potential at the origin. FigQ] shows the 
ratio between the frequencies of the second and the first 
mode UJ2/1-J1 as a function of 7, for the case of 20 ions. 
For N > 20 this ration does not depend on the number 
of ions neither on the value of A. For 7 = 1, and assum- 
ing small oscillations around the equilibrium positions, 
the ratio U2/0J1 is not well defined since the second mode 
does not exist because of the constant character of the 
force produced by the external potential. For 7 > 2 se- 
vere limitations in the storage capacity of the network, 
similar to those encountered in the harmonic case, ap- 
pear. The situation changes when 7 < 1. In this case 
the two lowest modes become nearly degenerated. Wc 
obtain larger number of stored patterns for 7 = 0.5 but 
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FIG. 3: Final overlap averaged over 500 initial configurations 
as a function of the initial overlap, for the two lowest vibra- 
tional modes of (a) 20 Ca + ions and (b) 40 Ca + ions in a po- 
tential V = A\x\ - 5 with A = 6.6 x W 20 J/m 1/2 . The black 
squares (triangles) correspond to the first (second) pattern. 



the desired degeneracy is only obtained for the first two 
collective modes. The results obtained with this optimal 
trapping potential are displayed in FigEl where the final 
overlap as a function of the initial overlap for the two 
lowest vibrational modes is depicted for (a) 20 Calcium 
ions and (b) 40 ions. The final overlap of the two lowest 
patterns and the corresponding ones with all the spins 
flipped is close to 1 up to two spin flips for N = 20. This 
shows that the system is able to recover the four pat- 
terns up to two spin flips. For N = 40 this holds up to 



eight spin flips. Note that the system does not always 
recover the initial pattern but instead sometimes reaches 
a slightly deformed configuration, which differs only in 
one spin flip from the original one. For this reason m/ 
in Fig. I^a) is not exactly 1 for one and two spin flips 
and up to eight in Fig. E^b). In Fig. E^a) the probabil- 
ity of recovering the initial pattern or one that differs in 
one spin flip with the initial one is 100% (96.9%) for the 
second (first) mode for the initial random spin flip of one 
spin and 98.7% (95.1%) for the second (first) mode for 
the initial random spin flip of two spins. In Fig. E^b) the 
probability of recovering of the two modes is above 98% 
up to three initial random spin flips and above 97% up 
to eight initial spin flips. 

Summarizing, we have shown that spin-ion systems can 
be used to implement NN models. We have calculated 
the storage capacity and robustness of such systems and 
showed how to control this storage capacity by changing 
the trapping potential and/or applying external longitu- 
dinal magnetic fields. Furthermore, such models should 
allow to implement quantum neural networks by apply- 
ing transverse magnetic fields that permit tunneling pro- 
cesses between stored patterns. Finally, NN provide a 
clear example of robust distributed information which 
could be exploited as a quantum memory. 
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